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INTRODUCTION 


The  Office  ol  Naval  Research  (ONR)  lias  requested  rhe  Naval  Civil 
Engineering  baboratory  (NCEL)  to  investigate  the  possibility  of  expe¬ 
diently  mod  fying  and  controlling  the  shape  of  the  beach.  As  part  of 
this  investigation,  NCEL  has  developed  a  numerical  model,  based  on  the 
physics  of  nearshore  sediment  transport  processes,  that  is  capable  of 
simulating  beach  profile  response,  to  waves. 

Quantitative  models  for  nearshore  sediment  transport  processes  are 
very  useful  for  predicting  nearshore  bathymetric  cuanges,  which  include 
changes  in  both  the  planform  shape  of  the  beach  and  the  beach  profile. 
Planform  changes  are  typically  associated  with  an  interruption  of  the 
longshore  littoral  drift  (e.g.,  by  a  groin).  Beach  profile  changes  are 
typically  associated  with  variations  in  the  cross-shore  sediment  trans¬ 
port  rate  which  is  driven  by  the  incident  wave  field.  Only  beach  pro¬ 
file.  changes  will  be  addressed  in  this  report. 

Beach  profile  changes  are  further  classified  according  to  time 
scale.  Short-term  events  range  from  a  few  hours,  when  associated  with  a 
single  storm  event,  to  many  months,  when  associated  with  seasonal  wave 
changes.  bong-term  events  typically  describe  a  rise  in  sea  level  and 
have  a  time  scale  measured  in  tens  of  years.  The  present  model  pertains 
only  to  short-term  events,  as  defined  above. 

EXISTING  SURFZONE  SEDIMENT  TRANSPORT  MODET.S 

Scientists  and  engineers  have  used  time-stepping  line  models  to 
simulate  the  bathymetric  evolution  of  the  coastal  zone  for  many  years 
(e.g.,  Pelnard-Considere,  1956).  These  models  represent  the.  nearshore 
bathymetry  by  one  or  more,  depth  contours.  These  contours  remain  constant 
as  the  model  evaluates  the  effects  of  waves  over  a  given  time  increment. 
At  the  end  of  the  time  increment,  bathymetric  changes  are  made.  The 
process  is  repeated  for  each  time  increment  for  all  sets  of  wave  condi¬ 
tions. 

The  first  to  include  cross-shore  sediment  transport  was  Bakker 
(1968),  who  developed  a  two- line  model  which  utilized  the  shoreline  and 
an  arbitrary  depth  contour  to  represent,  the  nearshore  bathymetry.  His 
transport  relationship  compares  the  cross-shore  separation  between  con¬ 
tour  lines  at  any  given  time  step  with  the  cross-shore  separation  at  the 
equilibrium  conditions  for  the  given  waves  The  offshore  transport,  rate 
per  unit  width,  Q^,  is: 

Qx  =  'ly[x,  -  (x„  -  w)  1  t  1 ) 
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where : 


~  rate  constant  with  dimensions  ( I,/T) 

X1’  X2  ~  cross-shore  positions  of  the  contours 

w  =  distance  (x^  -  x^)  for  the  equilibrium  profile 

Due  to  the  lack  of  a  transport  relationship  based  upon  the  physics 
of  sand  transport,  there  are  serious  shortcomings  in  this  approach. 

These  include:  (1)  a  lack  of  asymmetry  in  the  rate  of  offshore  and  on¬ 
shore  sediment  transport,  and  (2)  a  lack  of  accounting  for  wave  height, 
wave  period,  and  sediment  size  in  the  factors  w  and  q  . 

Swart  (1974)  used  the  general  form  of  Equation  1  in  bis  treatment 
of  offshore  transport  by  incorporating  both  sediment  and  wave  character¬ 
istics  in  his  constants.  Although  his  approach  to  both  the  equilibrium 
profile  and  the  proportionality  constant  were  more  sophisticated  than 
those  used  by  Bakker,  it  was  still  a  schematization  of  the  profile 
development  and  not  a  model  based  on  fundamental  sediment  transport 
mechan  ics . 

Stive  and  Battjes  (1984),  on  the  other  hand,  developed  an  offshore 
sediment  transport  model  based  on  an  equation  describing  the  time-averaged 
undertow  velocity.  It  predicts  the  offshore  transport  of  sediment  during 
storm  conditions  given  the  height  of  the  incident  waves  and  the  sediment 
load  in  the  near-bed  boundary  layer.  The  model  does  not  address  onshore 
transport . 

Recently,  many  cross -shore  sediment  transport  models  have  been 
developed  (Watanabo  and  Dibajnia,  1988;  Kuo  and  Chen,  1988;  Diegaard,  et 
al.,  1988;  and  Stive,  1988)  that  address  both  onshore  and  offshore  trans¬ 
port.  The  application  of  these  models,  though,  tends  to  be  rather  com¬ 
plex  . 

This  research  presents  a  beach  profile  evolution  model  which 
retains  much  of  the  simplicity  of  Bakker's  two-line  development,  while 
incorporating  several  significant  improvements.  Rather  than  using  a 
schematization  of  the  profile  development  to  simulate  the  cross-shore 
sediment  flux,  the  model  utilizes  cross-shore  sediment  transport  rela¬ 
tionships  (onshore  and  offshore)  based  on  the  fluid/sediment  mechanics 
in  the  surfzone.  Using  a  time-stepping  routine,  the  model  simulates  the 
nearshore  beach  profile  changes  by  applying  these  relationships. 

MODEL  DEVELOPMENT 

For  modeling  purposes,  the  beach  is  divided  into  two  regions,  the 
nearshore  region  which  extends  from  the  berm  crest  to  the  inflection 
point,  and  the  offshore  region  which  extends  from  the  inflection  point 
to  the  depth  of  closure  (see  Figure  1).  When  the  sediment  transport 
rate  is  negative,  the  beach  is  eroded  and  sediment  is  transferred  to  the 
offshore  region.  When  the  transport  rate  is  positive,  sand  is  removed 
from  the  offshore  region  and  deposited  in  the  nearshore  region. 

The  beach  profile  evolution  model  is  based  on  three  equations:  (1) 
the  continuity  oration,  (2)  a  cross-shore  sediment  transport  relationship 
developed  by  Bailard  (1981),  and  (3)  a  depth-of-cut  relationship  as  defined 
by  Hallermeier  (1977). 
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Conservation  of  Mass  Equation 


The  conservation  of  mass  requirement  states  that  sand  is  neither 
lost  nor  gained  by  the  control  volume  enclosing  the  Leach  profile.  This 
relationship  can  be  expressed  as: 


3  V 

at 


+ 


IS  = 

3x 


0 


(2) 


where:  V  =  total  volume  of  the  beach  profile  per  unit  width 

Q  =  transport  rate 
x  =  cross -shore  coordinate 


Cross-Shore  Sediment  Transport  Equation 


Baila’-d's  (1981)  cross-shore  sediment  transport  equation  is  based 
on  the  energetics  concepts  developed  by  Bagnold  (1963).  Bailard  gener¬ 
alized  Bagnold's  unidirectional  stream  transport  equation  to  allow  for 
time-varying  flow  over  an  arbitrary  sloping  bottom.  The  resulting 
equation  for  the  time-averaged  cross-shore  sediment  transport  rate, 

<  i  > ,  is : 
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where : 


density  of  water 

drag  coefficient  for  the  bed 

efficiency  factors 


tan  ,  W 


factors  based  on  sediment  characteristics 


tan  B  =  bed  slope 

Vr  * 

u^,  u  =  orbital  velocity  moments 

u  =  primary  component  of  the  oscillatory  velocity  field 
^1’  ^2  =  wave  skewness  parameters 


6  =  undertow 

u 
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Using  the  Nearshore  Sediment  Transport  Study  field  data  sets  (Gable, 

1979;  Gable,  i'<80),  spatially-averaged  empirical  estimators  were  devel¬ 
oped  for  the  wave  velocity  field,  its  moments,  the  skewness  parameters, 
and  the  undertow  (Bailard,  1985).  The  resulting  linear  equations  require 
only  the  incident,  significant  wave  height,  beach  slope,  and  median  sedi¬ 
ment  diameter  as  input. 

In  principle.  Equation  3  represents  a  balance  of  the  offshore- 
directed  transports  imposed  by  the  downslope  component  of  gravity  and 
the  time-averaged  undertow  with  the.  onshore  directed  transport  imposed 
by  the  asymmetry  of  the  wave  orbital  velocity  field.  At  equilibrium, 
these  transports  exactly  balance,  and  the  beach  is  in  dynamic  equilibrium 
with  the  incident  wave  conditions. 

Initial  Conditions 

Dean  and  Maurmeyer  (1983)  note  that  beach  response  time  is  as  much 
as  an  order  of  magnitude  larger  than  the  duration  of  most  hydrologic 
conditions.  Therefore,  conditions  seldom  persist  long  enough  to  com¬ 
pletely  dominate  the  development  of  a  beach  profile  bringing  it  to 
equilibrium.  More  typically,  a  beach  profile  is  the  result  of  a  long 
history  of  wave  conditions.  This  suggests  that  the  initial  conditions 
of  a  response  model  must  accurately  convey  the  true  state  of  the  beach 
to  successfully  simulate  its  subsequent  response. 

The.  user  is  required  to  input  the  following  values  to  define  the 
initial  conditions  of  the  beach  profile: 

•  Berm  height  -  a  constant  value  defined  as  the  vertical 
distance  from  the  still  water  line  to  the  crest  of  the  berm. 

•  Berm  slope  -  an  average  value  defining  the  ratio  of  the  height 
of  the  berm  to  the  horizontal  distance  between  the  mean 
shoreline  and  the  crest  of  the  berm. 

•  Foreshore  slope  -  the  initial  slope  of  the  submarine  beach 
profile  in  the  nearshore  zone. 

•  Offshore  slope  -  the  initial  slope  of  the  beach  profile  in  the 
offshore  zone,  excluding  the  horizontal  portion  of  the  bar. 

•  Depth  of  slope  break  -  the  initial  depth  from  the  still  water 
line  to  the  point  where  the  foreshore  slope  meets  the  offshore 
s  1  ope . 

•  bength  of  bar  -  the  initial  length  of  the  longshore  bar. 

Boundary  Conditions 

The  absolute  shoreward  boundary  of  the  model  is  the  termination  of 
the  negative  x-axiu.  This  boundary  is  far  shoreward  of  the  mean  shoreline 
which  acts  as  the  practical  shoreward  boundary  of  the  system.  The  mean 
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shoreline  has  one  degree  of  freedom  (x-axis),  and  its  movement  is  deter¬ 
mined  by  the  accretion  or  erosion  of  sand  in  the  -'earshore  region.  The 
berm  is  of  fixed  height  and  slope.  Its  position  is  uniquely  determined 
by  that  of  the  mean  shoreline. 

The  seaward  boundary  of  the  model  is  defined  by  a  fixed  depth  of 
closure.  This  is  the  depth  at  which  transport  ceases  for  the  median 
annual  significant  wave  (Hallermeier ,  1981).  That  is: 

D  =  1.418  (H  )(T  )»Ji7SD 

u  mm’ 

where:  DQ  =  depth  of  negligible  transport,  meters 

H  =  annual  median  wave  height,  meters 
m 

T  =  annual  median  wave  period 
m 

g  =  acceleration  of  gravity,  meters/second 
D  =  median  grain  size,  millimeters 

Note,  this  is  not  the  ultimate  depth  of  no  transport  for  all  wave  condi¬ 
tions,  but  merely  a  convenient  average  of  the  anticipated  values,  which 
will  vary  with  both  wave  height  and  wave  steepness.  The  user  is  required 
to  input  an  annual  median  deep  water  significant  wave  height  and  wave 
period  so  that  the  model  can  define  this  boundary. 

The  horizontal  position  of  the  seaward  boundary  is  determined  by 
the  location  at  which  the  depth  of  closure  falls  on  the  initial  offshore 
s  lope . 

Inflection  Point 

In  the  past,  modelers  (Bakker,  1968;  Swart,  1974)  have  used  a  line 
cf  constant  depth  to  define  the  division  between  the  nearshore  and  off¬ 
shore  regions.  In  the  present  model,  this  inflection  point  is  a  function 
of  both  the  incident  wave  height  and  the  time  history  of  the  profile. 

The  inflection  point  has  two  degrees  of  freedom  (x,z).  Its  movement  is 
defined  as  a  function  of  the  ratio  of  the  volume  of  sand  transported, 

Q  ,  to  the  total  volume  of  sand  transported  within  a  given  wave  regime, 

Qfc  (see  Figures  2  and  3).  The  incremental  movement  of  the  inflection 
point,  d,,  is  defined  as: 


d.  =  (D  -  D , )Q  /Q 
1  c  1  x  t 

where:  d.  =  incremental  movement 

1  at  time  step  "i" 

=  equilibrium  depth  of 

=  initial  depth  of  cut 


(4) 

of  the  inflection  point 

cut 

for  a  given  wave  regime 
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The  equilibrium  depth  of  cut  is  defined  by  Hallermeier  (1977)  as 
the  limiting  depth  of  erosion.  T t  is  a  functior  of  the  incident  wave 
height  and  wave  steepness: 

0  =  2.23  H  -  68 . 5( H  2/gT  2)  (5) 

c  s  s  s 

where:  H  =  significent  wave  height,  meters 

~  significant  wave  period 

The  horizonta1  position  of  the  inflection  point  is  determined  by 
the  location  at  which  the  depth  of  cut,  as  defined  by  Equation  5,  falls 
on  the  nearshore  slope. 


MODEL  APPLICATION 
Offshore  Transport 

The  offshore  transport  mode  occurs  when  the  present  nearshore  beach 
slope  is  steeper  than  the  equilibrium  slope  for  the  given  set  of  wave 
conditions.  In  this  mode,  sand  is  transported  from  the  nearshore  region, 
the  beach  slope  flattens,  and  the  berm/mean  shoreline  retreats.  Sand 
accumulates  in  the  offshore  region  forming  a  bar  whose  depth  is  determined 
by  the  inflection  point  and  whose  width  is  determined  by  the  volume  of 
sand  transported. 

Figure  4  shows  the  evolution  of  a  summer  beach  profile  that  is 
subjected  to  larger  waves.  A  wave  histogram  is  presented  in  Figure  5. 

The  steep,  reflective  nearshore  slope  flattens  as  the  berm/mean  shore¬ 
line  retreats  and  the  inflection  point  moves  offshore.  As  the  beach 
continues  to  reform  due  to  the  new  wave  conditions,  the  volume  transport 
rate  decreases  and  the  bar  grows  in  smaller  increments  as  it  moves  more 
slowly  offshore. 

Onshore  Transport 

The  onshoie  transport  mode  occurs  when  the  present  nearshore  beach 
slope  is  flatter  than  the  equilibrium  condition.  Sand  is  transported 
from  the  offshore  region  to  the  nearshore  as  the  nearshore  beach  steepens 
and  the  berm  accumulates  sediment. 

Figure  6  shows  the  evolution  of  a  winter  beach  profile  that  is  sub¬ 
jected  to  benign  conditions.  A  wave  histogram  is  presented  in  Figure  7. 
The  offshore  bar  begins  to  diminish  in  size  and  move  shoreward  as  the 
inflection  point  moves  to  shallower  water.  Sand  accumulates  on  the  berm 
forcing  the  mean  shoreline  to  advance.  A  steep,  reflective  nearshore 
slope  is  formed.  The  bar  continues  to  diminish  in  size  and  move  shore¬ 
ward,  much  like  a  ridge  and  runnel  system.  After  the  bar  is  welded  to 
the  teach  face,  the  nearshore  slope  continues  to  steepen  as  it  moves 
slowly  toward  the  equilibrium  condition. 


6 


Storm  Simulation 


Both  wave  da^a  and  beach  profile  data  were  collected  at  Leadbetter 
Beach  in  Santa  Barbara,  California  during  the  Nearshore  Sediment  Trans¬ 
port  Study.  This  data  set  includes  the  beach  profile  response  to  a 
storm.  Using  the  input  conditions  provided  by  the  prestorm  surveys  and 
the  wave  data  set  (.see  Figure  8),  it  was  possible  to  simulate  the  beach 
profile  response  using  the  cross-shore  sediment  transport  model. 

Figure  9  shows  the  measured  beach  profiles.  The  profiles,  from  the 
crest  of  the  berm  down  to  the  minus  2-meter  mark,  were  measured  via  wading 
profiles  throughout  the  storm  (Gable,  1980).  The  offshore  portion  of 
the  profile  was  measured  with  a  fathometer  before  the  storm  (Moore,  1982). 

Note  that  in  the  fieid  profiles  shown  in  Figure  9,  sand  was  removed 
from  the  beach  profile  and  not  just  redistributed  (as  would  be  in  keeping 
with  the  model  stipulation  of  conservation  of  mass).  Closer  inspection 
reveals  that  the  volume  of  sand  lost  came  from  the  nearshore  region 
(roughly  corresponding  to  the  surf  scone).  This  suggests  that  the  long¬ 
shore  transport  (driven  by  the  longshore  current  which  is  confined  to 
the  surfzone)  was  not  constant  in  space,  but  rather  has  resulted  in  a 
net  erosion  for  the  measured  beach. 

Since  this  beach  profile  response  model  does  not  take  into  account 
sediment  removed  by  longshore  transport,  compensation  must  be  made  for 
this  net  loss  of  sand  sc  the  field  profiles  can  be  compared  with  the 
model  simulation.  To  accomplish  this,  the  final  field  profile  was 
shifted  seaward  to  add  back  the  sand  that  was  lost  to  longshore  current- 
induced  erosion,  assuming  a  triangular  distribution  intersecting  the 
seafloor  at  the  point  of  closure  (s^e  Figure  10).  This  shifted  profile 
was  then  compared  to  that  predicted  by  the  model  (see  Figure  11). 

The  model  compares  reasonably  well  with  the  field  profile,  with  the 
exception  that  it  has  overpredicted  the  erosive  action  of  the  storm.  As 
shown  in  Figure  11,  the  length  of  the.  bar  and  the  depth  of  the  slope 
break  are  well  predicted,  yet  slightly  overpredicted.  The  flattening  of 
the  nearshore  beach  slope  and  the  retreat  of  the  mean  shoreline  are  over¬ 
predicted  due  to  the  excess  erosion. 

SUMMARY  AND  CONCLUSIONS 

The  present  model  shows  promise  as  a  method  for  simulating  cross¬ 
shore  sediment  transport.  It  is  a  simple  scheme  that  uses  only  significant 
wave  height,  wave  period,  sediment  size,  and  initial  beach  profile  as 
input.  Nevertheless,  the  mod-  1  effectively  simulates  erosion  as  well  as 
accretion.  The  model  tends  toward  an  equilibrium  profile  given  sufficient 
time  and  a  constant  wave  climate,  but  in  most  applications  this  is  never 
realized.  It  handles  events  that  have  time  scales  ranging  from  hours  to 
years . 

All  of  the  major  characteristics  of  beach  response  are  simulated. 
Offshore  transport  occurs  very  rapidly  during  storms  causing  development 
of  a  simple  barred  profile.  Onshore  movement,  on  the  other  hand,  is  a 
very  slow  process  in  relation  to  offshore  transport.  When  a  bar  is 
present  it  depicts,  in  a  rather  simple  fashion,  the  shoreward  marching 
of  a  ridge  and  runnel  system.  When  no  bar  is  present,  a  narrow,  steep 
beach  is  formed  in  response  to  the.  lower  waves. 
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Despite  its  ability  to  simulate  all  of  the  major  characteristics  of 
beach  response,  there  are  a  number  of  qua! ifications  concerning  this 
simple  model  that  must  be  emphasized  in  order  to  permit  its  useful 
application  as  a  coastal  engineering  tool: 

1.  The  model  assumes  a  conservation  of  mass.  This  implies  that 
the  longshore  component  of  sediment  transport  is  assumed  to  be  constant 
in  both  time  and  space.  This  is  often  not  the  case.  When  1ongshore 
transport  is  present,  it  may  be  appropriate  to  use  the  two-d Imens iona1 
cross-shore  transport  model  in  conjunction  with  a  two-dimensional  long¬ 
shore  sediment  transport  model  that  would  predict  the  sediment  gained  or 
lost  for  a  given  slice  cf  beach.  There  a^e  many  currently  available, 
for  example,  the  models  by  Le  M^Hauti  and  Soldate  (1980),  or  Dailard 
(1985)  would  be  adequate. 

2.  The  model  is  spatially-averaged  within  the  nearshore  region. 

This  means  that  the  cross-shore  distribution  of  the  sediment  transport 
rate  is  ignored.  This  eliminates  the  possibility  of  adequately  simulating 
multiple  bars. 

3.  The  offshore  region  acts  merely  as  a  source/sink,  dri/en  entirely 
by  the  dynamics  of  the  nearshore  zone. 

4.  The  empirically  derived  wave  velocity  moment  equations  are  limited 
to  the  relatively  narrow  range  of  conditions  found  on  the  southern  Californ 
coast  1 ine . 

5.  Due  to  a  bias  in  the  original  cross-shore  sediment  transport 
equa*  'on,  the  model  tends  toward  overpredicting  erosion  in  the  nearshore 
region  (as  defined  in  Figure  1). 
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1 1 


Figure  1.  Beach  schematic. 


of  the  inflection  point  during  erosion  of  the 
i  region. 


Figure  3.  Movement  of  the  inflection  point  during  accretion  in  the 
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Figure  4.  Evolution  of  the  beach  profile:  reflective  profile 
subjected  to  large  waves. 
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Figure  5.  Wave  histogram  for  evolution  of  the  beach  profile  due  to 
large  waves. 
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Figure  7.  Wave  histogram  for  evolution  of  the  beach  profile  due  to 
small  waves. 
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Figure  10.  Measured  beach  profiles,  Leadbetter  Beach,  CA.  Final 

profile  has  been  shifted  seaward  to  compensate  for  the  sand 
removed  by  the  longshore  transport. 
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Appendix  A 


COMPUTER  CODE  FOR  BFACH  PROFILE  RESPONSE  MODF.I. 


The  numerical  model,  CROSHORE,  is  a  two-dimensional  beach  profile 
response  model.  It  requires  the  user  to  supply  an  input  file,  entitled 
WAVES. DTA  (see  Appendix  B) ,  that  contains  wave  and  wave  condition 
duration  information.  All  other  input  information  that  is  required  is 
specifically  requested  by  CROSHORE  at  the  beginning  of  a  run.  All 
requested  information  must  be  supplied. 

The  numerical  model  "CROSHORE"  was  written  in  BASIC,  specifically, 
TurhoBasic  which  is  a  compiled  basic  program  published  by  Borland. 

Since  most  of  the  commands  in  TurhoBasic  are  similar  to  those  in 
interpretive  basic,  translation,  it  necessary,  should  be  fairly  easy. 

To  facilitate  the  person  who  wishes  to  "get  inside"  of  the  model,  every 
effort  was  made  to  add  comments  describing  each  variable.  Also,  comments 
describe  the  overall  function  of  each  loop. 

The  prc  gram  was  written  for  use  on  an  IBM-AT  compatible  personal 
computer.  All  of  the  applications  described  in  this  report  were  run  on 
a  Zenith  286  (AT  compatible.)  personal  computer  with  a  math  coprocessor. 
The  longest  application  took  less  than  2  minutes. 
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BEACH  PROFILE  RESPONSE  SIMULATION  MODEL  -  "CROSHORE" 


This  program  begins  with  an  initial  beach  profile  and 

wave  conditions  and  determines  the  beach  response  to  the  waves. 

All  units  used  in  this  program  are  metric  (meters -ki logram-second) . 

ESTABLISH  INITIAL  CONDITIONS 


"ENTER  MEDIAN  GRAIN  DIAMETER  IN  MILLIMETERS" 

D 

"ENTER  INITIAL  FORESHORE  SLOPE  IN  DEMICAL  FORM" 

Tanfo 

"ENTER  INITIAL  OFFSHORE  SLOPE  IN  DECIMAL  FORM" 

Tanof  f 

"ENTER  INITIAL  DEPTH  OF  SLOPE  BREAK" 

Dsb 

"ENTER  INITIAL  BAR  LENGTH" 

Lhar 

"ENTER  HEIGHT  OF  BERM" 

Dbf 

"ENTER  SLOPE  OF  BERM  FACE  IN  DECIMAL  FORM" 

Tanbf 

"ENTER  NUMBER  OF  WAVE  CONDITIONS  TO  BE  RUN" 

N 

"ENTER  ANNUAL  MEDIAN  WAVE  HEIGHT  AND  PERIOD" 

Hm ,  Tm 

” i F  OUTPUT  AT  CHANGE  OF  WAVE  CONDITIONS  ONLY,  ENTER  l" 

A 

•1  THEN  INPUT  "ENTER  NUMBER  OF  HOURS  BETWEEN  OUTPUTS" , Honrs 

'DRAG  COEFFICIENT 
'COEFFICIENT  OF  BEDLOAD  TRANSPORT 
'COEFFICIENT  OF  SUSPENDED  LOAD  TRANSPORT 
'SPECIFIC  GRAVITY  OF  SAND 
'ACCELERATION  OF  GRAVITY  (meters/sq-sec) 
'AT  REST  VOLUME  CONCENTRATION  OF  BEACH 
'ANGLE  OF  REPOSE  OF  SAND  IN  DECIMAL  FORM 
'TIME  STEP  INTERVAL  IN  SECONDS 


CLS 
PRINT 
IN  PIH’ 

PRINT 
INPUT 
PR  I  NT 
INPUT 
PRINT 
INPUT 
PRINT 
INPUT 
PRINT 
INPUT 
PR  I  NT 
INPUT 
PRINT 
INPUT 
PRINT 
INPUT 
PRINT 
INPUT 
IF  A< 

I 

Cf=0. 0009 

Eh~0 . 13 

F,s=0 . 032 

Ss=2.65 

G=9 .  R 

No-0 . 6 

Tana r=0 . 63 

Dt=3600 
1 

Quw-G*( Ss -  1 )*No 
W=(  8. 925/1000*I))*(SQR(  I+9S*(  Ss  -  1  )*D"  3)  -  1  ) 
Dr,=  1.41 8*Hmn*Tmn*SQR(  0/  (  5*1) ) ) 

Xc-I)sb/Tari  fo+(  Dc-Dsb)  /Tanof  f  +  Lbar 

t 


Xms 1 ~0 

Xbf-Xms I -Dbf /Tanbf 
Xsb=I)sb/Tan  fo 
KKK-1 


'SUBMERGED  UNIT  WEIGHT  OF  BEACH  SAND 
'FALL  VELOCITY  OF  SAND  GRAINS 
'DEPTH  TO  NEGLIGIBLE  TRANSPORT 
'DISTANCE  TO  Dc 


'DISTANCE  TO  THE  MEAN  SHORELINE 
'DISTANCE  TO  THE  CREST  OF  THE  BERM 
'DISTANCE  TO  THE  SLOPE  BREAK 
'GRAPHICS  COUNTER 
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WAVE  HEIGHT  AND  SPAN  OF  DURATION 


DIM  H(IOO),  S(IOO) 

OPEN  "WAVES. DTA"  FOR  INPUT  AS  #1 
OPEN  "PROFILE. DTA"  FOR  OUTPUT  AS  #3 
CALL  Grapher 
CALL  Plotter 


CALCULATIONS  LOOPS  -  TOTAL  VOLUME  TRANSPORTED  AND  BEACH  TRANSFORMATION 

t 

FOR  L=1  TO  N 

INPUT  in,  H(L),S(L) 

IF  A=1  THEN  Hours=S(L) 

Sumqt=0 

T  ime=0 

Dsb2=2*H(L) 

Knt=0 

1=0 
I 

'TOTAL  VOLUME  TRANSPORTED  LOOP 

I 

Qyt=i 

DO  UNTIL  Qyt<0 . 000 1 
Knt=Knt+l 

Lknt=0.001*Knt  'SMALLER  ITERATIVE  COUNTER 

B=0 . 1 -Lknt  'ITERATIVE  BEACH  SLOPE 

Psil=(2. 27*B-0. 141)*H(L)-I2. 7*B+0.547  ’WAVE  VELOCITY  MOMENT  1 

Psi2=(3.52*B-0.287)*H(L)-25*B+l. 13  'WAVE  VELOCITY  MOMENT  2 

Deltau=( -6 . 5*B-0 . 0158)*H( L) -4 . 29*B+0 .138  'MEAN  WATER  VELOCITY 

Um=( 28 . 5*B-0 . 288)*H( I»)  +  l . 29*B+0 .283  'PRIMARY  COMPONENT  OF  WAVE  VELOCITY 

Qybt=Cf*UrrT  3*Eb/(Quw*Tanar)*(Ps i 1+1 . 5*Del tau-0. 6*B/Tanar) 

Qyst=-Cf*Um"4*Es/(Quw*W)*(Psi2+2 . 4*Del tau-1 . 25*Es*Um*B/W) 

Qyt=Qybt+Qyst  'ADDITION  OF  BEDLOAD  AND  SUSUPENDED  LOAD  VOLUME  TRANSPORT  RATES 
LOOP 

Qt=ABS( (Dbf+0. 5*Dsb)*(Dsb/B-Dsb/Tanfo) )  'TOTAL  VOLUME  TRANSPORTED  FOR  EQUIL. 

t 

'BEACH  TRANSFORMATION  LOOP 

t 

T imestepp i ng: 

DO  UNTIL  ( I >=Hours  OR  Time>=S(L)) 

Sumq=0 

Psil=(2.27*Tanfo-0. 14 1 )*H(L) -  12 . 7*Tanfo+0. 547  'WAVE  VELOCITY  MOMENT  1 

Ps i 2=f 3 , S2*Tanfo-D. 287 )*JI( L) -25*Tanfo+l . 1 3  'WAVE  VELOCITY  MOMENT  2 

De 1 tau=( -6 . S*Tanfo-0 . 0 1 58)*H(  I.) -4 . 29*Tan fo+0 .138  'MEAN  WATER  VELOCITY 

Um=( 28. 5*Tanfo-0. 288)*H(L)+1 . 29*Tanof+0. 283  'PRIMARY  COMPONENT  OF  WAVE  VELOCITY 

Qyb=-Cf*Um  3*Eb/(Quw,vTanar  )*(  Ps  I  I  +  !  .  i;*Dp  I  tau-0 .  A,vTan  f o/Tanar ) 

Qsy=-C  f*Um  4*Es/(  Quw*W  )*(  Ps  i  2+2 . 4*De  1  tail  -  1 . 25*F,sUm*Tan  fo/W) 

Qy=+Qyb+Qys  'ADDITION  OF  BEDLOAD  AND  SUSPENDED  LOAD  VOLUME  TRANSPORT  RATER 

t 

Sumq=Sumq+Qy*I)t 


'SUM  OF  VOLUME  TRANSPORTED  FOR  H(L) 
'TIME  COUNTER 

'EQUILIBRIUM  SLOPE  BREAK  DEPTH  FOR  H(L) 
'TOTAL  VOLUME  TRANSPORTED  LOOP  COUNTER 
'BEACH  TRANSFORMATION  LOOP  COUNTER 


A-  3 


Sumqt=Sumqt+Qy*Dt 
1  =  1  +  1 

T  ime=T  icne+Dt  /  3600 


IF 


Lbar>0  OR  Sumqt>0  THEN 

Deldsb=(Dsb2-Dsb)*ABS( Sumq/Qt ) 

Dsb=Dsb+De ldsb 

Xsb=Xsb+Deldsb/Tanfo 

Des lxms l=-( Sumq-0 . 5*Deldsb*( 2*Xbar 

Xms  l=Xins  l+De  lxms  1 


POSITIVE  TRANSPORT  IS  OFFSHORE 
'CHANGE  IN  Dsb 
'NEW  Dsb 
'NEW  Xsb 

Deldsb*(Tanfo-Tanoff)))/(Dbf+0. 5*Dsb) 
'CHANGE  IN  Xms 1 
'NEW  Xms 1 


Xbf=Xms 1 -Dbf /Tanb  f 

De 1 lbar=( Deldsb/Tanof  f-De ldsb /Tan fo) +2*Sumq/( Dc -Dsb) 

Lbar=Lbar+Del lbar 

Tanfo=Dsb/(Xsb-Xmsl) 

Tanof f=(Dc-Dsb)/(Xc(Xsb+Lbar) ) 


'NEW  Xbf 
'CHANGE  IN  Lbar 
'NEW  Lbar 
'NEW  Tanfo 
'NEW  Tanoff 


IF  Lbar<0  THEN  Lbar=0 


'DEPTH  OF  PROFILE  INTERSECTION 

'CHANGE  IN  Xms 1 

'NEW  Xsml 

'CHANGE  IN  Xsb 

'NEW  Xsb 

'NEW  Xbf 

'NEW  Tanfo 

'NEW  Tanoff 


1=0 

CALL  Plotter 

IF  Time^SCL)  THEN  GO  TO  Timestepping 

NEXT  L 

CLOSE  #1 

CLOSE  # 3 

END 

t 

'SUB  PROGRAM  THAT  DRAWS  THE  GRAPHICS 

I 

SUB  Grapher 
screen  9 

SHARED  Xc , Dc , Db  f , Sw 1 , ENDX ,  ENDY 


ELSE 

Adi=l/2 

Bdi=-(Dbf+Dsb+Dc) 

Cdi=Dsb*Dc 

Di=( -Bdi-SQR(Bdi" 2-4*Adi*Cdi) )/( 2*Adi) 
Del xms l=-Sumq/(Dbf+Di/2) 

Xms l=Xms l+De lxms 1 
Delxsb=Sumq/(Dc-Di) 

Xsb=Xsb+De lxsb 
Xb f=Xms 1 -Dbf/Tanbf 
Tanfo=Dsb/ (Xsb=Xsml) 

Tanof  f=(Dc-Dsb)/(Xc-Xsb) 

END  IF 
LOOP 

I 
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DRAW  THE  BOX 


line  (80,56)-(600,  306), ,B 
» 

'PRINT  THE  TITLES 

t 

locate  24,30 

print  "CROSS  SHORE  DISTANCE  (m) 
LOCATE  1,31 

PRINT  "BEACH  PROFILE  CHANGES" 
LOCATE  3,11 

PRINT  "DEPTH  FROM  SWL  (m) 

f 

'X  AXIS  TICK  MARKS  AND  LABELS 

t 

ENDX=INT(Xc/ 100+1 )+2 
X=l 

DO  UNTIL  X=ENDX+ 1 
Xx=520*(X/F.NDX)+80 
LINE  (Xx,292)-(Xx,287) 
Labelx=(X-3)*100 
Xl=2*INT(X/2) 

IF  X1<>X  THEN 

LOCATE  22,65*(X-1) /ENDX+9 

PRINT  Labe lx 

END  IF 

X=X+1 

LOOP 


Y  AXIS  TICK  MARKS  AND  LABELS 


'MARKS  THE  POSITION  OF  THE  SWL 
'SAVES  THE  POSITION  OF  THE  SWL 


ENDY=INT( (Dbf )+(Dc+2) ) 

Y=1 

DO  UNTIL  Y=ENDY 

Yy=250*(Y/ENDY)+42 

Dcr=INT(Dc+2) 

Ly=(Dcr-ENDY) - 1 
Labely=-(Ly+Y) 

LINE  (80,Yy)-(85,Yy) 

IF  Labely=0  THEN 

LINE  ( 80 , Yy )  -  (600,Yy) 
Swl=Yy 
END  TF 

Y1=2*INT( Y/2) 
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IF  Y1<>Y  THEN 

LOCATE  1 7 . 5*(  Y- 1 )  /F.NDY+5 , 8 

PRINT  Labely 

END  IF 

Y=Y+1 

LOOP 

END  SUB 

t 

'SUB  PROGRAM  THAT  PLOTS  THE  PROFILES 

t 

SUB  Plotter 
screen  9 

SHARED  Xb f ,  Db f ,  Xms  1 ,  Sw  1 ,  Xsb ,  Dsb ,  Xc ,  Dc ,  Lbar ,  ENDX ,  F.NDY ,  KKK 

i 

'CONVERTING  TO  GRAPHICS  UNITS 
1 

Dbfp=Swl-250*(Dbf /ENDY) 

Dsbp=Swl+250*(Dsb/ENDY) 

Dcp=Sw 1+25  0* ( Dc/ENDY ) 

Xbfp=520*( (Xbf / 100+2 ) /ENDX+80 

Xms  lp=520*( ( Xms 1 / 100+2 ) /ENDX+80 

Xsbp=520*( ( Xsb/ 1 00+2 ) /ENDX+80 

Xbarp=520*( ( (Xsb+Lbar)/ 100+2 )/ENDX)+80 

Xcp=520*(  (Xc/  100+2) /ENI)X)+80 
» 

' PLOTTING 

f 

LINE  (Xbfp.Dbfp)  -  (Xms lp, Swl ) , KKK 
LINE  (Xmslp,Swl)  -  (Xsbp,Dsbp) ,KKK 
IF  Lbar>0  THEN 

LINE  (Xsbp,Dsbp)  -  (Xbarp.Dsbp) ,KKK 
LINE  (Xbarp,Dsbpl  -  (Xcp, Dcp) , KKK 

ELSE 

LINE  (Xsbp.Dsbp)  -  (Xcp, Dcp) , KKK 
END  IF 

IF  KKK=1  OR  KKK=1 3  THEN  Mask=&HAAA 
IF  KKK=5  THEN  Mask=&HBBB 
IF  KKK=9  THEN  Mask=&HCCC 
KKK=KKK+4 
END  SUB 
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INFORMATION  CONCERNING  INPOT  FILES 


INPUT  FILES 

Input  files  can  be  created  with  any  line  editor.  The  file, 
entitled  WAVES. DTA,  must  be  in  ASCII  format  with  two,  comma  delimited 
fields  per  line.  The  first  field  should  be  wave  height,  in  meters;  the 
second  field  should  be  duration  of  wave  conditions,  in  hours. 


WAVE MAKER  PROGRAM 

For  those  not  familiar  with  line  editors  and/or  file  structure,  a 
simple,  program,  WAVEMAKER,  has  been  written.  It  provides  a  simple,  but 
slightly  cumbersome,  method  of  creating  a  suitable  input  file,  WAVES. DTA. 
Note  that  each  time  WAVEMAKER  is  run  it  automatically  deletes  the 
previous  input  file.  To  avoid  deleting  previous  input  files  one  must 
rename  the  files  not  being  used  in  the  present  application. 
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’ "WAVEMAKER" 

t 

'This  program  inputs  and  stores  (and  edits,  minimally)  the  wave 

'conditions  required  to  run  the  program  "CROSHORE." 

! 

'JUST  IN  CASE  LOOP 

t 

CLS 

PRINT  "THIS  PROGRAM  WILL  ERASE  FILESPECS  WAVES. DTA  AND  WAVES. XXX" 
PRINT  "IF  THIS  IS  NOT  ALRIGHT  WITH  YOU,  ENTER  9" 

INPUT  Ers 
IF  F,rs=9  THEN 
GOTO  Bye 

ELSE 

KILL  "WAVES.*” 

END  IF 

t 

' INPUT  INSTRUCTIONS  LOOP 

I 

CLS 

PRINT  "ENTER  THE  NUMBER  OF  WAVE  CONDITIONS  TO  BE  USED" 

INPUT  N 

DIM  H(100),  S(100),  E(100) 

'DATA  INPUT 

t 

OPEN  "WAVES. DTA"  FOR  OUTPUT  AS  #1 
FOR  1=1  TO  N 

PRINT  "ENTER  WAVE  HEIGHT  AND  DURATION  OF  CONDITION"l  "OF"  N 
INPUT  H( I ) , S( I ) 

WRITE  #1 , H( I ) , S( I ) 

NEXT  I 
CLOSE  in 

I 

'CHECKING  DATA  BY  REVIEWING  SCREEN  OUTPUT 

t 

Pagain: 

PRINT  "PLEASE  CHECK  DATA  SET  FOR  ERRORS" 

PRINT  "  " 

PRINT  "CONDITION  HEIGHT  DURATION" 

OPEN  "WAVES. DTA"  FOR  INPUT  AS  #1 
FOR  J=t  TO  N 

INPUT  ifl,  H(J),S(J) 

PRINT  J,H( J) ,S( J) 

NEXT  J 
CLOSE  in 

PRINT  "IF  DATA  REQUIRES  EDITING,  ENTER  l" 

INPUT  Checkdat 
IF  Checkdat=l  THEN 
CALL  Ed i tor 
GOTO  Pagain 
END  IF 
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Bye : 

PRINT  "YOUR  DATA  SET  IS  NOW  STORED  IN  WAVES. DTA" 
END 


'SUB-PROGRAM  THAT  ALLOWS  EDITING  OF  DATA 
» 

SUB  Editor 
CLS 

SHARED  N,H(),S() 

PRINT  "DATA  WILL  BE  CHANGED  IN  ASCENDING  ORDER  OF  CONDITION  NUMBER" 

PRINT  "  " 

PRINT  "ENTER  NUMBER  OF  CONDITIONS  TO  BE  EDITED" 

INPUT  Ne 

NAME  "WAVES. DTA"  AS  "WAVES. XXX"  'PUTS  PREVIOUS  DATA  SET  INTO  A  DUMMY  FILE 
OPEN  "WAVES. XXX"  FOR  INPUT  AS  #2 
OPEN  "WAVES. DTA"  FOR  OUTPUT  AS  #3 
FOR  Cor.=  l  TO  Ne 

PRINT  "ENTER"  Con  "OF"  Ne  "CONDITIONS  TO  BE  EDITED 
INPUT  E(Con) 

NEXT  Con 
Knt=l 

PRINT  Knt , E(Knt) 

DELAY  10 

FOR  C=1  TO  N  'LOOP  WHERE  THE  EDITING  TAKES  PLACE 

INPUT  #2,  H(C) ,  S(C) 

IF  C=E(Knt)  THEN 
Knt=Knt+l 

PRINT  "WAVE  HEIGHT  =  "  H(C),  "DURATION  ="  S(C) 

PRINT  "ENTER  APPROPRIATE  DATA" 

INPUT  Hn(C) , Sn(C) 

H(C)=Hn(C) 

S(C)=Sn(C) 

END  IF 

WRITE  #3,  H(C) , S(C)  'PUTS  NEW  INFO.  INTO  OLD  FILE 

NEXT  C 
CLOSE  # 2 

CLOSE  #3  'ELIMINATES  THE  DUMMY  FILE 

KILL  "WAVES. XXX" 

END  SUB 
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DISTRIBUTION  LIST 


AF  /  6550  ABG/DER,  Patrick  AFB,  FL 

AF  /  CSR  4200  (K.  Davidson),  Patrick  AFB,  FL 

AFESC  /  TIC  Lib,  Tyndall  AFB,  FL 

ARMY  /  HQDA  (DAEN-ZCM) ,  Washington,  DC 

ARMY  /  POJED-O,  APO  San  Francisco, 

ARMY  BELVOIR  R&D  CEN  /  STRBE-AALO,  Ft.  Belvoir,  VA 

ARMY  BELVOIR  R&D  CEN  /  STRBE-FS,  Ft.  Belvoir,  VA 

ARMY  BELVOIR  R&D  CEN  /  STRBE-JB,  Ft.  Belvoir,  VA 

ARMY  ENGRG  DIST  /  Lib,  Seattle,  WA 

ARMY  ENGRG  DIST  /  Lib,  Portland,  OR 

ARMY  ENGRG  DIST  /  Phila,  Lib,  Philadelphia,  PA 

ARMY  EWES  /  Lib,  Vicksburg,  MS 

ARMY  EWES  /  WESCD  (TW  Richardson),  Vicksburg,  MS 
ARMY  EWES  /  WESCP-D  (Vallianos),  Vicksburg,  MS 
ATLANTIC  RICHFIELD  CO  /  RE  Smith,  Dallas,  TX 
BATTELLE  /  D.  Frink,  Columbus,  OH 
BATTELLE  /  New  Eng  Marine  Rsch  Lab,  Duxbury,  MA 
BECHTEL  CIVIL,  INC  /  K.  Mark,  San  Francisco,  CA 
BLAYLOCK  WILLIS  &  ASSOC  /  T  Spencer,  San  Diego,  CA 
BROWN  &  ROOT  /  Ward,  Houston,  TX 
BUREAU  OF  RECLAMATION  /  Smoak,  Denver,  CO 
CAL  STATE  UNIV  /  C.V.  Chelapati,  Long  Beach,  CA 
CALIF  MARITIME  ACADEMY  /  Lib,  San  Ramon,  CA 
CATALINA  MARINE  SCI  CEN  /  Lib,  Avalon,  CA 
CATHOLIC  UNIV  /  CE  Dept  (Kim)  Washington,  DC 
CHEVRON  OIL  FLD  RSCH  CO  /  Strickland,  La  Habra,  CA 
CHILDS  ENGRG  CORP  /  K.M.  Childs,  Jr.,  Medfield,  MA 
CINCLANTFLT  /  CE  Supp  Plans  Offr,  Norfolk,  VA 
CINCUSNAVEUR  /  London,  UK,  FPO  New  York, 

CLEVELAND  STATE  UNIV  /  CE  Dept  ( Akinmusuru) ,  Cleveland,  OH 

CNO  /  DCNO,  Logs,  OP-424C,  Washington,  DC 

CNO  /  DCNO,  Logs,  OP -45 2,  Washington,  DC 

COLLINS  ENGRG,  INC  /  M  Garlich,  Chicago,  IL 

COLORADO  STATE  UNIV  /  CE  Dept  (Criswell),  Ft.  Collins,  CO 

COLORADO  STATE  UNIV  /  CE  Dept  (W  Charlie),  Fort  Collins,  CO 

COLUMBIA  GULF  TRANSMISSION  CO  /  Engrg  Lib,  Houston,  TX 

COM  GEN  FMF  /  LANT,  SCE,  Norfolk,  VA 

COM  GEN  FMF  /  PAC,  FEO,  Camp  HM  Smith,  HI 

COMDT  COGUARD  /  G-ECV,  Washington,  DC  , 

COMFAIR  /  MED,  SCE,  Naples,  Italy,  FPO  New  York, 

COMFLEACT  /  SCE,  FPO  Seattle, 

COMNAVACT  /  PWO,  London,  UK,  FPO  New  York, 

COMNAVSURF  /  Code  N42A,  Norfolk,  VA 

COMNAVSURF  /  LANT,  CO,  Norfolk,  VA 

COMNAVSURF  /  Pac,  Code  N-4,  San  Diego,  CA 

CONSTRUCTION  TECH  LABS,  INC  /  G.  Corley,  Skokie,  IL 

CONTINENTAL  OIL  CO  /  0.  Maxson,  Ponca  City,  OK 

CORNELL  UNIV  /  Lib,  Ithaca,  NY 

DAMES  &  MOORE  /  Lib,  Los  Angeles,  CA 

DILLINGHAM  CONSTR  CORP  /  (HD&C),  F  McHale,  Honolulu,  HI 


DTRCEN  /  Code  1175,  Annapolis,  MD 

DTRCEN  /  Code  119,  Annapolis,  MD 

DTRCEN  /  Code  125,  Annapolis,  MD 

DTRCEN  /  Code  1568,  Annapolis,  MD 

DTRCEN  /  Code  4120,  Annapolis,  MD 

DUKE  UNIV  /  CE  Dept  (Muga),  Durham,  NC 

EASTPORT  INTL,  INC  /  JH  Osborn,  Mgr,  Ventura,  CA 

EDWARD  K  NODA  fit  ASSOC  /  Honolulu,  HI 

EFA-SW  /  Code  101.1,  San  Diego,  CA 

EG7G  WASH  ANAL.  SVC  CEN  /  Bonde,  Gaithersburg,  MD 

EWI  ENGINEERING  ASSOCIATES  /  Jack  Cox,  Middleton,  WI 

EXXON  PRODUCTION  RSCH  CO  /  Offshore  Op  Div,  Houston,  TX 

FLORIDA  ATLANTIC  UNIV  /  Ocean  Engrg  Dept  (Martin),  Boca  Raton,  FL 

FLORIDA  ATLANTIC  UNIV  /  Ocean  Engrg  Dept  (McAllister),  Boca  Raton,  FL 

FLORIDA  ATLANTIC  UNIV  /  Ocean  Engrg  Dept  (Su),  Boca  Raton,  FL 

FLORIDA  ATLANTIC  UNIV  /  Ocean  Engrg  Dept  (Wolf)  Boca  Raton,  FL 

FLORIDA  INST  OF  TECH  /  CE  Dept  (Kalajian),  Melbourne,  FL 

GEORGE  WASHINGTON  UNIV  /  Engrg  &  App  Sci  Schl  (Fox),  Washington,  DC 

GEORGIA  INST  OF  TECH  /  CE  Schl  (Mazanti),  Atlanta,  GA 

GLIDDEN  CO  /  Rsch  Lib,  Strongsville,  OH 

INST  OF  MARINE  SCIENCES  /  Dir,  Morehead  City,  NC 

INST  OF  MARINE  SCIENCES  /  Lib,  Port  Aransas,  TX 

INTL  MARITIME,  INC  /  D.  Walsh,  San  Pedro,  CA 

IOWA  STATE  UNIV  /  CE  DEPT,  Ames,  IA 

KLIEGER,  PAUL  /  CE,  Northbrook,  IL 

LAMONT- DOHERTY  GEOL  OBSERVATORY  /  McCoy,  Palisades,  NY 

LEHIGH  UNIV  /  CE  Dept,  Hydraulics  Lab,  Bethlehem,  PA 

LEHIGH  UNIV  /  Marine  Geotech  Lab,  Bethlehem,  PA 

LONG  BEACH  PORT  /  Engrg  Dir  (Allen),  Long  Beach,  CA 

MAINE  MARITIME  ACADEMY  /  Lib,  Castine,  ME 

MARATHON  OIL  CO  /  Gamble,  Houston,  TX 

MARCORBASE  /  Base  Maint  Dept,  Camp  Lejeune,  NC 

MARCORBASE  /  Maint  Offr,  Camp  Pendleton,  CA 

MARITECH  ENGRG  /  Donoghue,  Austin,  TX 

MARITIME  ADMIN  /  MMA,  Lib,  Kings  Point,  NY 

MC  CLELLAND  ENGRS ,  INC  /  Lib,  Houston,  TX 

MICHIGAN  TECH  UNIV  /  CO  Dept  (Haas),  Houghton,  MI 

MIT  /  R.V.  Whitman  Cambridge,  MA 

MOBIL  R&D  Corp  /  Offshore  Engrg  Lib,  Dallas,  TX 

MOFFATT  &  NICHOL  ENGRS  /  R.  Palmer,  Long  Beach,  CA 

MUESER  RUTHLEDGE,  CONSULTING  ENGRS  /  Lib,  New  York,  NY 

NAS  /  NI,  Code  183,  San  Diego,  CA 

NATL  ACADEMY  OF  ENGRY  /  Alexandria,  VA 

NAVACAD  /  Kriebel,  Annapolis,  MD 

NAVCOASTSYSCEN  /  CO,  Panama  City,  FL 

NAVEODTECHCEN  /  Tech  Lib,  Indian  Head,  MD 

NAVFACENGCOM  /  Code  04A1D,  Alexandria,  VA 

NAVFACENGCOM  /  Code  04A3,  Alexandria,  VA 

NAVFACENGCOM  /  Code  04A3C,  Alexandria,  VA 

NAVFACENGCOM  /  Code  04A4E,  Alexandria,  VA 

NAVFACENGCOM  /  Code  06,  Alexandria,  VA 

NAVFACENGCOM  /  Code  07,  Alexandria,  VA 

NAVFACENGCOM  /  Code  163,  Alexandria,  VA 


NAVFACENGCOM  /  Code  183  (Gebhart),  Alexandria,  VA 
NAVFACENGCOM  CHESDIV  /  Code  FP0-1C,  Washington,  DC 
NAVFACENGCOM  LANTDIV  /  BR  OFC,  DIR,  Naples,  Italy,  FPO  New  York, 
NAVFACENGCOM  LANTDIV  /  Code  411,  Norfolk,  VA 
NAVFACENGCOM  PACDIV  /  Code  2011,  Pearl  Harbor,  HI 
NAVFACENGCOM  SOUTHDIV  /  Code  04A3,  Charleston,  SC 
NAVFACENGCOM  SOUTHDIV  /  Code  102B,  Charleston,  SC 
NAVFACENGCOM  WESTDIV  /  Code  09P/20,  San  Bruno,  CA 
NAVFACENGCOM  WESTDIV  /  Code  102,  San  Bruno,  CA 
NAVFACENGCOM  WESTDIV  /  Code  1833,  San  Bruno,  CA 
NAVFACENGCOM  WESTDIV  /  Code  406.2  (Smith),  San  Bruno,  CA 
NAVFACENGCOM  WESTDIV  /  Code  407,  San  Bruno,  CA 

NAVFACENGCOM  WESTDIV  /  PAC  NW  Br  Offc,  Code  C/42,  Silverdale,  WA 

NAVOCEANO  /  Code  6200  (M  Paige),  NSTL,  MS 

NAVOCEANO  /  Lib,  NSTL,  MS 

NAVOCEANSYSCEN  /  Code  94,  San  Diego,  CA 

NAVPGSCOL  /  E.  Thornton,  Monterey,  CA 

NAVPHIBASE  /  SCE,  San  Diego,  CA 

NAVSEASYSCOM  /  Code  05M3,  Washington,  DC 

NAVSHIPREFAC  /  Lib,  FPO  San  Francisco, 

NAVSHIPYD  /  Tech  Lib,  Portsmouth,  NH 
NAVSTA  /  CO,  Long  Beach,  CA 
NAVSTA  /  Engr  Div,  PWD,  FPO  Miami, 

NAVSTA  /  PWO,  Mayport,  FL 

NAVSTA  PANAMA  CANAL  /  Code  54,  FPO  Miami, 

NAVSUPPACT  /  PWO,  Holy  Loch,  UK,  FPO  New  York, 

NAVWPNSTA  EARLE  /  Code  092,  Colts  Neck,  NJ 
NBS  /  Bldg  Mat  Div,  Mathey,  Gaithersburg,  MD 
NETPMSA  /  Tech  Lib,  Pensacola,  FL 

NEW  YORK  STATE  MARITIME  COLLEGE  /  Longobardi,  Bronx,  NY 

NIEDORODA,  AW  /  Gainesville,  FL 

NMCB  /  3,  Ops  Offr,  FPO  San  Francisco, 

NMCB  /  40,  CO,  FPO  San  Francisco, 

NOAA  /  Joseph  Vadus,  Rockville,  MD 

NOARL  /  Code  440,  NSTL,  MS 

NORDA  /  Code  112 ISP,  NSTL,  MS 

NORDA  /  Code  352,  NSTL,  MS 

NORDA  /  Code  440,  NSTL,  MS 

NORTHWEST  ENGRG  CO  /  Grimm,  Bellevue,  WA 

NORTHWESTERN  UNIV  /  CE  Dept  (Dowding),  Evanston,  IL 

NRL  /  Code  6127,  Washington,  DC 

NSC  /  SCE,  Norfolk,  VA 

NUSC  DET  /  AUTEC  Andros  Is,  Asst  Wpns  Offr,  FPO  Miami, 

NUSC  DET  /  Code  2143  (Varley),  New  London,  CT 
NUSC  DET  /  Code  44  (Munn),  New  London,  CT 
NUSC  DET  /  Code  TA131,  New  London,  CT 
NUSC  DET  /  PWO,  New  London,  CT 
NWSTAT  /  R0ICC ,  Colts  Neck,  NJ 
NY  CITY  COMMUNITY  COLLEGE  /  Lib,  Brooklyn,  NY 
OCNR  /  Code  1121  (EA  Silva),  Arlington,  VA 
0ICC  /  Engr  and  Const  Dept,  APO  New  York, 

OREGON  STATE  UNIV  /  CE  Dept  (Bell),  Corvallis,  OR 
OREGON  STATE  UNIV  /  CE  Dept  (Yim),  Corvallis,  OR 


OREGON  STATE  UNIV  /  Oceanography  Scol,  Corvallis,  OR 

PACIFIC  MARINE  TECH  /  M.  Wagner,  Duvall,  WA 

PENNSYLVANIA  STATE  UNIV  /  Applied  Rsch  Lab,  State  College,  PA 

PHIBCB  /  1,  CO,  San  Diego,  CA 

PILE  BUCK,  INC  /  Smoot,  Jupiter,  FL 

PURDUE  UNIV  /  Engrg  Lib,  West  Lafayette,  IN 

PWC  /  CO,  Oakland,  CA 

PWC  /  Code  123C,  San  Diego,  CA 

SAN  DIEGO  STATE  UNIV  /  CE  Dept  (Noorany),  San  Diego,  CA 

SEATECH  CORP  /  Peroni,  Miami,  FL 

SHANNON  &  WILSON,  INC  /  Lib,  Seattle,  WA 

SHELL  OIL  CO  /  Boaz,  Houston,  TX 

SHELL  OIL  CO  /  E.  Doyle,  Houston,  TX 

SUBASE  /  SCE,  Pearl  Harbor,  HI 

SUPSHIP  /  Tech  Lib,  Newport  News,  VA 

TAMPA  PORT  AUTHORITY  /  Engrg  Dept  (Schrader),  Tampa,  FL 
TECHNOLOGY  UTILIZATION  /  K  Willinger,  Washington,  DC 
TEXAS  A&M  UNIV  /  CE  Dept  (Herbich),  College  Station,  TX 
TEXAS  A&M  UNIV  /  Ocean  Engr  Pro j ,  College  Station,  TX 
UCT  /  TVO,  CO,  Port  Hueneme,  CA 
UCT  ONE  /  CO,  Norfolk,  VA 

UNIV  OF  ALASKA  /  Biomed  &  Marine  Sci  Lib,  Fairbanks,  AK 

UNIV  OF  CALIFORNIA  /  CE  Dept  (Fourney),  Los  Angeles,  CA 

UNIV  OF  CALIFORNIA  /  CE  Dept  (Gerwick),  Berkeley,  CA 

UNIV  OF  CALIFORNIA  /  CE  Dept  (Mitchell),  Berkeley,  CA 

UNIV  OF  CONNECTICUT  /  Lib,  Groton,  CT 

UNIV  OF  DELAWARE  /  CE  Dept,  Newark,  DE 

UNIV  OF  HAWAII  /  Manoa,  Lib,  Honolulu,  HI 

UNIV  OF  HAWAII  /  Ocean  Engrg  Dept  (Ertekin) ,  Honolulu,  HI 

UNIV  OF  ILLINOIS  /  Metz  Ref  Rm,  Urbana,  IL 

UNIV  OF  MICHIGAN  /  CE  Dept  (Richart),  Ann  Arbor,  MI 

UNIV  OF  TENNESSEE  /  CE  Dept  (Kane),  Knoxville,  TN 

UNIV  OF  WASHINGTON  /  RL  Terrel,  Edmonds,  WA 

UNIV  OF  WISCONSIN  /  Great  Lakes  Studies  Cen,  Milwaukee,  WI 

US  GEOLOGICAL  SURVEY  /  Marine  Geological  Offc,  Reston,  VA 

USCINCPAC  /  Code  J44,  Camp  HM  Smith,  HI 

USNA  /  Ocean  Engrg  Dept,  Annapolis,  MD 

VALLEY  FORGE  CORPORATE  CENTER  /  Franklin  Research  Center,  Norristown,  PA 

VIATEUR  DE  CHAMPLAIN  /  Inst  of  Maritime  Engrg,  Matane,  Quebec 

VIRGINIA  INST  OF  MARINE  SCI  /  Lib,  Gloucester  Point,  VA 

VSE  /  Ocean  Engrg  Gp  (Murton),  Alexandria,  VA 

WOODWARD-CLYDE  CONSULTANTS  /  R.  Cross,  Oakland,  CA 

WOODWARD -CLYDE  CONSULTANTS  /  West  Reg,  Lib,  Oakland,  CA 

WORCESTER  POLYTECH  INST  /  J.M.  Sullivan,  Worcester,  MA 
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Air  pollution 

OCEAN  ENGINEERING 
Seafloor  soils  and  foundations 

Seafloor  construction  systems  ana  operations  (including 
diver  and  manipu'ator  tools) 
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